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ABSTRACT 

In the high density cores of globular clusters, multibody interactions are expected to be common, with the 
result that black holes in binaries are hardened by interactions. It was shown by Sigurdsson & Hernquist (1993) 
and others that 10 M Q black holes interacting exclusively by three-body encounters do not merge in the clusters 
themselves, because recoil kicks the binaries out of the clusters before the binaries are tight enough to merge. Here 
we consider a new mechanism, involving four-body encounters. Numerical simulations by a number of authors 
suggest that roughly 20-50% of binary -binary encounters will eject one star but leave behind a stable hierarchical 
triple. If the orbital plane of the inner binary is strongly tilted with respect to the orbital plane of the outer 
object, a secular Kozai resonance, first investigated in the context of asteroids in the Solar System, can increase 
the eccentricity of the inner body significantly. We show that in a substantial fraction of cases the eccentricity 
is driven to a high enough value that the inner binary will merge by gravitational radiation, without a strong 
accompanying kick. Thus the merged object remains in the cluster; depending on the binary fraction of black 
holes and the inclination distribution of newly-formed hierarchical triples, this mechanism may allow massive 
black holes to accumulate through successive mergers in the cores of globular clusters. It may also increase the 
likelihood that stellar-mass black holes in globular clusters will be detectable by their gravitational radiation. 

Subject headings: black hole physics — (Galaxy:) globular clusters: general — gravitational waves — stellar 
dynamics 



1. INTRODUCTION 

Globular clusters are outstanding testbeds for dynamics. As 
dense systems with ages many times their core relaxation time, 
they display such features as core collapse and mass segrega- 
tion, and they are almost certainly affected strongly by the pres- 
ence of even a small number of binaries. It has long been spec- 
ulated that various processes might produce relatively massive 
black holes in their cores (e.g., Wyller 1970; Bahcall & Os- 
triker 1975; Frank & Rees 1976; Lightman & Shapiro 1977; 
Marchant & Shapiro 1980; Quinlan & Shapiro 1987; Porte- 
gies Zwart et al. 1999; Ebisuzaki et al. 2001). Recent obser- 
vations of some dense clusters provide tentative evidence for 
black holes as massive as 2500 Mq at their centers (Gebhardt 
et al. 2000). 

Qualitatively, it seems entirely reasonable that large black 
holes should grow in the cores of many clusters. Even at birth, 
black holes are much more massive than the average star in a 
cluster, and hence they sink rapidly towards the core. When in 
the core, they tend to exchange into binaries. If the binary is 
hard (i.e., if its binding energy exceeds the average kinetic en- 
ergy of a field star), then a subsequent interaction with a field 
star tends to harden the binary (e.g., Heggie 1975). If this pro- 
cess is repeated often enough, the binary becomes tight enough 
that it can merge by gravitational radiation, and the black hole 
becomes larger. If these binaries merge while still in the cluster, 
sources in globulars could be excellent prospects for detection 
by the upcoming generation of gravitational wave instruments. 

However, it has been shown (e.g., Sigurdsson & Hernquist 
1993; Portegies Zwart & McMillan 2000) that if all black holes 
have initial masses of 10 Mq, three-body encounters alone do 
not lead to the formation of a large black hole at the center. 
The reason is that hardening in a binary-single interaction is 



accompanied by recoil, which kicks the binary out of the clus- 
ter before it can merge. Without additional effects, this means 
that the mergers occur well away from their host globulars. If 
the initial mass of a black hole is > 50 Mq, as may result from 
a high-mass low-metallicity star or rapid merger of main se- 
quence stars, it has enough inertia to remain in the core and 
grow by coalescence (Miller & Hamilton 2002). But what if 
only low-mass black holes are produced? 

Here we propose a new mechanism for the coalescence of 
low-mass black holes in globular clusters, involving binary- 
binary interactions. Studies of such four-body encounters have 
been comparatively rare, but have shown that in roughly 20- 
50% of the interactions the final state is an unbound single star 
plus a stable hierarchical triple system (Mikkola 1984; McMil- 
lan, Hut, & Makino 1991; Rasio, McMillan, & Hut 1995). This 
allows an important new effect: studies of planetary and stel- 
lar systems have shown that if there is a large relative incli- 
nation between the orbital planes of the inner binary and the 
outer object of the triple, then over many orbital periods the rel- 
ative inclination periodically trades off with the eccentricity of 
the inner binary, sometimes leading to very high eccentricities 
(Kozai 1962; Harrington 1968, 1974; Lidov & Ziglin 1976). 
In turn, this can enhance the gravitational radiation rate enor- 
mously, leading to merger without a strong kick and allowing 
even low-mass binary black holes in globulars to be potential 
gravitational wave sources. 

In § 2 we discuss the principles of this resonance, as de- 
rived in the case of three objects of arbitrary mass by Lidov 
& Ziglin (1976). To their treatment we add, in § 3, a simple 
term that accounts for general relativistic pericenter precession. 
We show that although, as expected, this precession decreases 
the maximum attainable eccentricity for a given set of initial 
conditions, the decrease is typically minor and thus there is sig- 
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nificant phase space in which the eccentricity resonance leads 
to rapid merger. In § 4 we use these results in a simple model 
for the mergers of black holes, and show that, depending on the 
fraction of black holes in binaries, this effect can lead to a dra- 
matic increase in the retention of black holes in globulars, and 
to the growth of ~ 10 2-3 M black holes in their cores. 



2. PRINCIPLES OF THE KOZAI RESONANCE 

When looking for changes in the orbital properties of a three- 
body system that extend over many orbital periods of both the 
inner binary and the outer tertiary, it is convenient to average the 
motion over both these periods, a procedure called double av- 
eraging. A general analysis of the double-averaged three-body 
problem has been performed to quadrupolar order for Newto- 
nian gravity by Lidov & Ziglin (1976) in Hill's case, in which 
the distance of the outer object (of mass m 2 ) from the inner 
binary (with component masses mo and mi < m ) is much 
greater than the semimajor axis of the inner binary. They find 
that for any set of three masses there is always a relative incli- 
nation of orbits such that an inner binary with arbitrarily small 
initial eccentricity will evolve to e = 1. For example, in the 
restricted three-body problem in which m 3> m 2 S> mi (e.g., 
the Sun, Jupiter, and an asteroid interior to Jupiter's orbit; see 
Kozai 1962), a relative orbital inclination of 90° will cause the 
asteroid to evolve to e = 1 in a finite time. 

However, the growth to such high eccentricity depends on 
a long series of perturbations from the tertiary that add co- 
herently, and hence requires certain phase relations. An extra 
source of precession of the pericenter can interfere with this. 
For example, the orbits of the moons of Uranus are tipped by 
97° with respect to Uranus' orbit around the Sun, but their ec- 
centricities stay relatively low due to precession introduced by 
the quadrupole moment of Uranus. In the case of black holes 
or other close massive objects, a similar role may be played by 
the effects of general relativity, which to lowest order includes 
precession of the pericenter. How does this affect the maximum 
eccentricity for a given set of initial conditions? 

Hill's approximation allows us to treat the system as two 
nested binaries: the inner pair composed of m and mi, and 
a second pair consisting of i) an object of mass m + mi lo- 
cated at their center of mass and ii) m 2 . Defining variables 
as in Lidov & Ziglin (1976), we let Mi = m + mi and 
M 2 = m + mi + m 2 be the total masses of the two binaries, 
and iii = m mi/Mi, and fi 2 = m 2 Mi/M 2 be their reduced 
masses. Let the semimajor axes and eccentricities of the two 
binaries be a\, e\ and a 2 , e 2 , and define i\ and i 2 to be the in- 
clinations of the binaries relative to the invariant plane of angu- 
lar momentum of the system. Finally, let \i — GmQm\m,ilM\ 
and e = 1 — e\. 

The double-averaged Hamiltonian H admits several integrals 
each of which yields a constant of the motion. First, the double- 
averaging procedure guarantees that ai and a 2 are constant. We 
keep terms in the Hamiltonian up to linear order in a\ /a 2 ; these 
quadrupolar terms dominate the evolution of the system for the 
high relative inclinations of interest here (see Ford, Kozinsky, 
& Rasio 2000). To this order, e 2 is also constant. The prob- 
lem has two constants of the motion that are related to angular 
momentum: 



e 1 / 2 cosii + /3cosz 2 




The constant (3 (a combination of the constants ai, a 2 , and 
e 2 ) represents the total angular momentum of the outer binary 
while a is the total system angular momentum (with contribu- 
tions from both the inner and outer binaries). Both a and (i are 
made dimensionless by dividing by L\ = ni^/GMiai, the an- 
gular momentum that the inner binary would have if it were on 
a circular orbit. 

The Hamiltonian, Ti, itself is constant. For convenience, we 



define U = ~k(W 
and obtain: 
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-2e + e cos 2 i" + 5(1 - e) sin 2 oj (cos 2 / - l) 



(2) 



which is equation (30) from Lidov & Ziglin (1976). Here u) is 
the argument of pericenter of the inner binary and the scaled 
angular momenta a, j3, and y/e form a triangle from which the 
relative inclination I = i 1 + i 2 can be obtained using the law 
of cosines: 

a 2 -/? 2 -6 

cos/= Wi ' (3) 

The maximum e (and hence minimum ei) occurs for oj = 0, and 
the minimum e (and hence maximum ei) occurs for oj = ir/2; 
see Lidov & Ziglin (1976). Therefore, given initial values for eo 
and ujq, the maximum eccentricity may be derived from conser- 
vation of W at lu = 7r/2. The time required to push the system 
from its minimum to maximum eccentricity, is of order 
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(e.g., Innanen et al. 1997), where 6 2 = a 2 -\/l — e\ is the 
semiminor axis of the tertiary and typically / <~ few for / near 
90°, which is the case of interest here. 

3. THE KOZAI RESONANCE WITH GR PRECESSION 

Post-Newtonian precession may be included in a couple of 
equivalent ways. One is to modify the Hamiltonian directly, 
by changing the gravitational potential to simulate some of the 
effects of general relativity. The modification of the potential 
is not unique, and depends on which aspect of general relativ- 
ity is to be reproduced (see Artemova, Bjornsson, & Novikov 
1996). For our purposes it is the precession of pericenter that 
is important (as opposed to, e.g., the location of the innermost 
stable circular orbit), and hence the correct lowest-order modi- 
fication is -GM/r -> (-GM/r)(l + ZGM/rc 2 ) (Artemova 
et al. 1996). Averaging the correction term over the orbits of the 
tertiary and inner binary, we obtain a correction to the double- 
averaged Hamiltonian of 
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This result may also be obtained from the first-order general 
relativistic precession rate of doj — (QttGMi/ [ai(l — e 2 )c 2 ] 
over one binary period (see Misner, Thorne, & Wheeler 
1973, p. 1110) using the equation of motion dco/dt = 
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-(2ky/e/Lt)(dW/dz) derived by Lidov & Ziglin (1976). Sub- 
stituting and integrating, we find that the first-order post- 
Newtonian contribution to W is 
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in agreement with equation (|5j). We have also checked our ex- 
pressions with direct numerical three-body integrations; note 
that equation corrects a factor of two error in equation (19) 
of Lin etal. (2000). 

Adding the new term Wpn to equation (Q) and making use 
of equation (||), we find 

W = -2e + ecos 2 / + 5(l-e) sin 2 w(cos 2 / - 1) +6» PN /e 1/2 . 

(7) 

As in the previous section, for a given set of initial conditions, 
one can therefore solve for the minimum e (maximum e), by 
setting u) = ir/2 and using the conservation of W. In general 
we expect that initially the inner binary will have significant 
eccentricity caused by perturbations during the four-body en- 
counter, but for simplicity we will assume that the initial ec- 
centricity is small enough that eo ~ 1 . In the restricted three- 
body problem in which mo S> m 2 mi and the initial rel- 
ative inclination is Iq, the approximate solution for e m ; n when 
5 cos 2 Jo « 3 (high inclination) and 6>pn <C 3 (weak preces- 
sion) is 
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iq ^> 6*pn this reduces to the Newtonian 
solution, in which the maximum eccentricity is e max = 
y/l - (5/3) cos 2 7 (Innanen et al. 1997). If instead I » tt/2 
so that 60 cos 2 I -C #pn> then e max ~ 1 — 9p N /9. More gen- 
erally, for any set of masses, if e — > 1 is allowed in the New- 
tonian problem then e max = 1 — 0(9p N ) when general rela- 
tivistic precession is included. Numerically, for Mi = Mq and 
ai =1 AU, 6» PN = 8 x 10- 8 (A/i/m 2 )(6 2 /ai) 3 . Equation (|) 
shows that in the restricted three-body problem the maximum 
possible eccentricity (minimum e m ; n ) is attained for the initial 
condition Jo = 7r/2 (initially perpendicular circular orbits). If 
mi has non-negligible mass, so that m 2 dominates the total 
angular momentum less, then the critical Iq increases (Lidov 
& Ziglin 1976). Figure 1 shows the critical inclination in the 
Newtonian case (#pn = 0) for several mass ratios and semima- 
jor axes. 

We want to know whether this process can cause the inner 
binary to reach a high enough eccentricity that it merges by 
gravitational radiation before the next encounter with a star in 
the globular cluster (which will typically alter the eccentricities 
and inclinations significantly). Encounters with black holes in 
globular clusters are usually dominated by gravitational focus- 
ing instead of the pure geometrical cross section; this is true 
within ~ 100AU of a 10 Mq black hole, where we have as- 
sumed a velocity dispersion of 10 km s _1 for the interlopers 
(see Miller & Hamilton 2002). In this limit the encounter time 
is 

r cnc » 6 x 10 B n^(l AU/a 2 ) (10 M /M 2 ) yr , (9) 

where the number density of stars in the core of the globular is 
10 6 n6 pc~ 3 . Note that it is the semimajor axis of the outermost 
object, m 2 , that sets the encounter time scale, because in a sta- 
ble hierarchical triple a 2 must be a factor of several greater than 



The timescale for merger by gravitational radiation for a high 
eccentricity orbit is (Peter 1964) 
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The steep dependence on eccentricity means that shrinkage of 
the orbit is dominated by the time spent near maximum eccen- 
tricity. Assuming that tgr S> t ovo i so that orbital decay occurs 
over many Kozai oscillation cycles, one finds that the fraction 

1/2 

of time spent near e max is of order (Innanen et al. 1997, 
equation [5]), so that r GR « 5 x 10 17 (^) (^j) 4 e 3 nlin yr. 
The condition for merger before an encounter is then simply 

r GR < T cnc- 
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Critical relative inclinations for evolution from e ss 
in the Newtonian case; attaining e = 1 requires a = (3 in equation 
the inner binary is composed of equal mass stars, mo = mi, and the labels on 
the curves indicate the mass ratio m2/m\. As the fraction of the total angu- 
lar momentum supplied by the tertiary increases (larger j3 and therefore larger 
"12 /mi or 02/01), the critical inclination trends toward 90°. 

Note that in the Newtonian case 6p^ = 0, all systems 
with the same masses, & 2 /ai, and Iq are dynamically iden- 
tical, in that the maximum eccentricity does not depend on 
the individual values of £> 2 and ai. The introduction of post- 
Newtonian precession breaks this scaling. If 6 2 /ai is fixed, 
then #PN a aj" 1 and therefore the maximum eccentricity at- 
tained is given by e m i n oc a^ 2 (cf. equation ^ for the restricted 
problem). The merger time is tgr, oc afe^ in oc a^ 2 . That is, 
a wider binary can be pushed to higher eccentricities, and actu- 
ally merge faster, than a closer binary. Note, however, that the 
solid angle for this orientation is proportional to #pn oc aj" 1 , 
because the optimum angle is usually close tow/2, so the solid 
angle is proportional to the cosine of the inclination. There- 
fore, if binary-binary interactions leave the binary and tertiary 
inclinations randomly oriented with respect to each other then 
a smaller fraction of wide binaries will fall into the optimal ori- 
entation. Qualitatively this means that as the binary is hard- 
ened by various interactions, every time a triple is formed it has 
a chance to push the eccentricity high enough that the binary 
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merges before the next encounter. The smaller the system, the 
larger the probability of such an orientation, because both the 
solid angle and the encounter time are larger. 

One way to quantify the probability of merger through the 
increase of eccentricity is to plot, as a function of the semima- 
jor axis of the inner binary, the range of relative inclinations 
such that merger occurs before the next encounter of a field 
black hole with the tertiary (which, being on a wide orbit, will 
interact before the inner binary will on average). In Figure 2, 
we assume three 10 M@ black holes, with a given a\ and a^. 
From <22 and an assumed number density of stars in the cluster 
(n = 10 6 pc -3 ), we compute the average time T cnc to the next 
encounter within a distance 02 of the system. We then deter- 
mine the range of initial inclinations / such that tqr < r enc , 
by solving for e m j n using equation ([7]) with the initial conditions 
ei = e2 = 0.01 and u> — 0. Note that for wider tertiary orbits, 
the total angular momentum of the system is dominated more 
by the tertiary (larger /3), and hence the relative inclination that 
gives the smallest possible £,„;„ is closer to 90° (see Figure 1). 
If a single Kozai oscillation cycle is longer than r cnc the system 
never attains the required high eccentricity. This causes the cut- 
off in the 02 = 10ai and a-i = 20cii curves in Figure 2; similar 
cutoffs exist at ai > 10 AU for the remaining two curves. 

1 I — 1 — 1 — 1 — 1 — I — 1 — r~n — 1 — gr~ 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 
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FlG. 2. — Inclination ranges for merger by gravitational radiation. For this 
graph, we assume that all three black holes have mass M = 10 Mq, and we 
assume a globular core number density n = 10 6 pc — 3 for calculating r cnc . 
The shaded regions indicate ranges of the relative inclination for which merg- 
ers will occur for each of four value of the semimajor axis ratio 02/01 : 3, 5, 10, 
and 20. The peaks, which occur at the locations predicted by the bottom curve 
in Fig. 1 , are truncated where the time to increase the eccentricity of the inner 
binary is greater than the mean time to an encounter (t cvo i < r cnc ). For com- 
parison, the dotted lines are the boundaries of the regions if general relativistic 
precession is suppressed; only for 02 = 20ai is there a noticeable difference. 



4. CONCLUSIONS 

The level of importance of the Kozai mechanism depends on 
several factors including: i) details of the interactions between 
two binaries, ii) details of the interactions between a triple, and 
either a binary or a single star, and iii) the fraction of black 
holes in binaries, which in turn relies on the iv) dynamics of the 
cluster itself. Understanding these interactions statistically will 
require extensive long-term simulations. However, the Kozai 
mechanism has the potential to be the dominant process in the 
interactions of stellar-mass black holes in globulars, if most 
such black holes are in binaries. When only three-body inter- 
actions are considered, very few black holes are retained by 
the clusters (only 8% in the simulations of Portegies Zwart & 
McMillan 2000). This occurs because the same processes that 
harden a binary toward an eventual merger also impart velocity 
kicks on the binary that ultimately eject it from the globular be- 
fore it can merge. In contrast, the majority of black holes can 
be retained if binary-binary interactions dominate. 

For example, suppose that a third of those interactions pro- 
duce stable triples. Subsequent interactions of the tertiary with 
field stars will change its eccentricity and semimajor axis. If 
the pericenter distance of the tertiary is less than a few times 
a\, then the triple system becomes unstable, normally by eject- 
ing its least massive member. Suppose that there are typically 
^2 encounters before the triple is disrupted in this way, and 
that each encounter of the tertiary that does not create an un- 
stable triple produces a new relative inclination I that is drawn 
from a uniform distribution in cos I. Suppose also that every 
time the inner binary interacts strongly its semimajor axis is 
decreased by ~ 20% (typical for strong interactions of three 
equal-mass objects; see, e.g., Heggie 1975; Sigurdsson & Phin- 
ney 1993). Then, in an n = 10 6 pc ~ 3 cluster there is a w 50% 
chance that the inner binary will merge before it hardens to 
ai ~ 0.2 AU, at which point the binary recoil velocity w rC coii 
exceeds the ^50 km s _1 escape speed typical of the cores of 
globulars (Webbink 1985). In an n — 10 5 pc ~ 3 cluster, en- 
counters are less frequent and the fraction rises to w 70%. 

Thus, depending on the binary fraction and other properties 
of black holes in globulars, the majority of black holes could 
merge before being ejected, and growth of intermediate-mass 
black holes in globulars may proceed naturally even if no black 
hole is formed with M > 10 M Q . This could influence stellar 
dynamics in the core, and the gravitational wave signals from 
globulars, and should be included in future simulations. 

This work was supported in part by NASA grant NAG 5- 
9756 and by NSF grant 5-23467. 
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